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Abstract 


A spatial template is important to study nearby supernova remnants (SNRs). For SNR G332.5-5.6, we report a 
Gaussian disk with a radius of about 1°06 to be a potential good spatial model in the y-ray band. Employing this 
new Gaussian disk, its GeV lightcurve shows a significant variability of about seven sigma. The y-ray observations 
of this SNR could be explained well either by a leptonic model or a hadronic model, in which a flat spectrum for 


the ejected electrons/protons is required. 
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1. Introduction 


The energy released by supernova events is of great 
importance for our understanding of the physics of the 
interstellar medium (ISM). Depending upon the environments, 
the supernova remnants (SNRs) can display a vast range of 
shapes (Whiteoak & Green 1996). The released energy comes 
into the electrons and protons, which interact with the interstellar 
radiation or protons, and finally emit the observable electro- 
magnetic signals up to high-energy bands (>100 MeV) (Acero 
et al. 2016). 

In the Anglo Australian Observatory (AAO)/United King- 
dom Schmidt Telescope (UKST) Ha survey, a new Galactic 
SNR was uncovered (Parker et al. 2005). This source with an 
unusual morphology, dubbed the “paperclip,” was discovered 
from the original Ho survey films, which was further confirmed 
as SNR G332.5-5.6 (Stupar et al. 2007). SNR G332.5-5.6 
shows three patches of filamentary emission with a total size of 
about 30' at the radio band (Reynoso & Green 2007). It has an 
extended X-rays morphology in the center region, which has a 
good correlation with radio emission detected at different 
frequencies (Suárez et al. 2015). SNR G332.5-5.6 is found to 
be a distance about 3.4 kpc with an age of 7—9 kyr (Zhu et al. 
2015). For its y-ray observations, such as in the Fermi High- 
Latitude Extended Sources (FHES) catalog, SNR G332.5-5.6 is 
the potential association of FHES J1642.1-5428, which has a 
uniform disk with a radius of about 0757 (Ackermann et al. 
2018). In an incremental version of the fourth full catalog of 
Fermi-LAT sources (4FGL-DR3), SNR G332.5-5.6 associates 
with an extended source of 4FGL J1642.1-5428e, which shares 
a disk radius of about 0°70 (Abdollahi et al. 2022). 

Morphology in the y-ray band is very different from that in 
the low-energy bands, such as radio or X-ray bands. Since the 
FHES catalog employed only 8 yr of observations by Fermi 


Large Area Telescope (Fermi/LAT), the spatial model of SNR 
G332.5-5.6 could be studied in deep employing more gamma- 
ray observations. We thus could perform morphology study in 
the high-energy gamma-ray bands, such as 14.3 yr Fermi-LAT 
observations. The paper is organized as follows. In Section 2, 
data analysis is presented. In Section 3, the physical origin of 
GeV emission of SNR G332.5-5.6 is explored. A summary and 
conclusions are presented in Section 4. 


2. Fermi/LAT Data Analysis 
2.1. Data Selection 


FermiPy version 1.2 package ' is used to analyze the Fermi / 
LAT data, which is a python package that facilitates analysis of 
the Fermi-LAT data with the latest Fermi Science Tools 
version 2.2.0? (Wood et al. 2017). Pass 8 events are selected in 
the region of interest (ROI) of 15? around the position of SNR 
G332.5-5.6, that is R.A., decl. = 2507538, —54°477 (Ajello 
et al. 2017). We chose the photons in the energy range between 
100 MeV and 1 TeV within 2008 August and 2022 December. 
The events were selected with event class of 128 and event type 
of 3. Events with zenith angles larger than 90? are excluded to 
avoid contamination from the Earth's limb. We use the 
standard data quality selection criteria (DATA. QUAL > 0) 
&&(LAT CONFIG == 1). We accepted a pixel size of 
0°02 spatially and eight logarithmic energybins per decade 
when performing the energy selection. PSR3. SOURCE V3 
instrument response function is used in our data analysis. 


! https: //fermipy.readthedocs.io/en/v1.2/ 


https: //github.com /fermi-lat /Fermitools-conda /wiki/Installation- 
Instructions 
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Table 1 
Spatial Analysis Results between 100 MeV and 1 TeV 
Spatial Template? R.A. Decl. Reg Ak —log(Likelihood) AAIC 
(deg) (deg) (deg) 

NONE Zr Ue 0 —16292446.11 0 

PS 250.925 + 0.063 —54.388 + 0.065 4 —16292463.69 —27.16 
GLEAM Zr EE e 2 —16292481.79 —67.36 
DISK 250.564 + 0.065 —54.245 + 0.075 0.536 + 0.037 5 —16292557.78 —213.34 
GAUSS 250.173 + 0.131 —54.401 + 0.117 1.061 + 0.153 5 —16292571.08 —239.94 


Note. 


è Spatial template for SNR G332.5-5.6, NONE represents null source, GLEAM is an extended map by radio observation, GAUSS is the Gaussian disk, PS is a single 


point source, DISK is the uniform disk. 


2.2. Basic Model 


A basic model is built by adding the background gamma- 
ray sources. 4FGL-DR3 is employed, in which we include all 
sources that 25? within SNR G332.5-5.6 center, as well as the 
diffuse galactic interstellar emission (gll iem, v07.fits) and 
the isotropic emission (iso PSR3 SOURCE V3 VI.txt). The 
spectral parameters of the diffuse galactic interstellar emis- 
sion, the isotropic emission and all sources that 3° within 
SNR G332.5-5.6 center are allowed to be free. We found there 
are 11 point sources as shown in Figure 1, whose spectral 
parameters are allowed to be free. Lastly, we removed the 
catalog source of 4FGL J1642.1-5428e in the ROI center and 
then built a basic model (hereafter the model with the spatial 
template of NONE). 


2.3. Spatial Analysis 


We evaluate the significance of all test spatial models, which 
is quantified through the likelihood ratio test with two statistical 
methods. One is the test statistic (TS) to compute the 
significance for a test model over the NONE spatial model, 
which is defined as 


TS = 2 x (In Lest — In Lyone), (1) 


Where Lies, and Lyone are the maximum likelihood values of 
the test model and of the basic NONE model respectively 
(Mattox et al. 1996). The likelihood value depends on the 
model and the observations. For the Fermi-LAT observations, 
the photon numbers in the ROI are 9065267, 1691719, 56799, 
2052 in the energy bands of 0.1-1 GeV, 1-10GeV, 
10-100 GeV and 100-1000GeV respectively. The other 
one is the Akaike information criterion (AIC) test to find 
out a better model between two nest models or unnest 
models (Akaike 1974; Tibaldo et al. 2018; Tang et al. 2021; 
Guo & Xin 2024), which is defined as 


AIC = 2k — 2inL (2) 
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Figure 1. TS map (7? x 7?) of the NONE model with pixel size of 0?1 between 


0.1 and 1000 GeV centering at SNR G332.5-5.6. Sources marked by white 
cross are 11 background point sources. 


where k is the number of degrees of freedom (d.o.f) and £ is the 
maximum likelihood value of a model. Usually, a model with a 
smaller AIC is the better one among several models. In this 
case, we employed the difference in AIC (AAIC) between two 
models to claim the prefer model, which could be calculated by 


AAIC = 2Ak — 2Aln Z. (3) 


Four spatial templates are tested with a single power-law 
(PL) spectral model, as shown in Table 1. TS maps excluding 
SNR G332.5-5.6 in each model is plotted in Figure 2. 
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Figure 2. TS maps (5° x 5°) after removing SNR G332.5-5.6 in four spatial templates. Top left: DISK template with a radius of ~0°54. Top right: GAUSS template 
with a radius of ~1°06. Bottom left: PS template. Bottom right: GLEAM template. The green circle represents the corresponding radius, green cross is the point source 


and green contour is from the radio observations. 


1. PS model. It shares the spatial template of a point source, scanned the disk radius from 0°01 to 1°20 using the tool of 


which has four additional free parameters (two for the 
free position and two PL spectral parameters). The 
resultant maximum TS is 35.16 and the best-fit position is 
about R.A., decl. = 250°925, —54°388. 

. GLEAM model. It is derived from radio observations of 
Galactic and Extragalactic Allsky MWA Survey between 
170 and 231 MHz is employed (Wayth et al. ), see 
Figure A1, which has two additional free parameters (two 
PL spectral parameters). The resultant maximum TS 
is 71.36. 


source finding. We found a best-fit radius of 07536 with 
the maximum TS of 223.34. The best-fit centered position 
is about R.A., decl. = 2507564, —54°245. 


. GAUSS model. It is a two-dimensional Gaussian disk, 


which also has five additional free parameters. The scanning 
method is same as that in the DISK model. The best-fit 6896 
containment radius is found at about 17061 with the 
maximum TS of 249.94. 'The best-fit centered position is 
about R.A., decl. = 2507173, —54°401. 


Results of AIC test are shown in Table 1. As seen, the 
GAUSS model is confirmed to be the best representation of the 
data, with a minimum AAIC of —239.94 with respect to the 


3. DISK model. It is a uniform disk, which has five additional 
free parameters (One for free radius, two for the free 
centered position and two PL spectral parameters). We 
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Figure 3. Lightcurve of SNR G332.5-5.6 with 14 time bins in the energy band 
of 0.1—1000 GeV. Data points are plotted with errors at confidence levels 95%. 
The light-blue shadows are the TS values. The red line is the constant flux 
employing the full time observation of 14.3 yr. 


NONE model. Considering that GAUSS also has the largest TS 
value among four test spatial models, the GAUSS model is the 
preferred model in our spatial analysis. In the latest Fermi-LAT 
catalog (4FGL-DR4), there are three SNRs with spatial 
template of GAUSS, such as SNR G51.3+0.1, SNR G150.3 
+4.5 and SNR G2922-0.5? Therefore GAUSS spatial 
template is selected for SNR G332.5-5.6 to produce its 
lightcurve (LC) and the spectral energy distribution (SED) in 
the following analyses. 


2.4. Temporal Analysis 


To test whether GeV emission of SNR G332.5-5.6 is 
variable temporally, we construct its lightcurve in GAUSS 
spatial model with a PL spectral model. We separate Fermi- 
LAT observations into 14 time intervals in the energy range 
from 0.1 GeV to 1 TeV. Energy flux in each time interval is 
calculated, which is plotted in Figure 3. The variability index 
(TS,ar) is defined as 

N 

TS = -27 In Li(Feonst) (4) 

i Li (Fi) 
where £; is the value of the likelihood corresponding to the i-th 
bin, F; is the best-fit flux for bin i, and Feonst is the best-fit flux 
for the full time assuming a constant flux (Nolan et al. 2012; 
Mo et al. 2023). The resultant TS,,. is 83.40 with 6.9c for the 
LC of 14 time bins, which suggests that there has a significant 
variability for the photon fluxes of SNR G332.5-5.6. 


2.5. Spectral Analysis 


Gamma-ray emission of the analyzed source is represented by a 
PL spectral model in the 4AFGL-DR3 catalog. To fully investigate 
the y-ray properties of this source, we considered other types of 
spectral models, such as broken PL (BPL), log-parabola (LP), and 
PLSuperExpCutoff (PLEC), which however do not result in any 


5. https://fermi.gsfc.nasa.gov/ssc/data/access/lat/ 14yr. catalog/ 
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Table 2 
Best-fit Physical Parameters in Leptonic Model and Hadronic Model 
Physical Model logio(A) Q/ Qj Ecutoff w° 
(GeV!) (TeV) (1058 erg) 
IC 41.06 2.01 1.05 3.08 
PD+Secondary IC 52.88 2.03 29.82 132.46 


Note. 
^ The total energy of particles (electrons or protons) in units of 107? erg. 


significant likelihood of improvement. We thus employed a PL 
spectral model to fit the GeV emission of SNR G332.5-5.6, e.g., 
<x El, The resultant photon index (Ly) is 2.14 + 0.04. The 
spectral energy distribution is plotted in Figure 4. 


3. Physical Model 


We explore the physical origins utilizing Fermi/LAT 
observations of SNR G332.5-5.6. Naima package? is 
employed to find the best-fit parameters in two physical 
models, such as a leptonic model and a hadronic model 
(Zabalza 2015). 


3.1. Leptonic Model 


We considered a leptonic model to fit observed GeV SED of 
SNR G332.5-5.6, e.g., the y-ray emission is generated through 
inverse Compton (IC) scattering of soft photons by relativistic 
electrons (Kamae et al. 2006; Meyer et al. 2010; Tang 2018). 
The ejected electronic distribution is assumed to follow an 
exponential-cutoff power-law function (PLEC), 


E\* E 
PM i A(Z) is (- | n 


where A is the normalization, Ep is the pivot energy fixed at 
1 TeV, a, is the electron injection spectral index, and Eeutoff 
stands for the cutoff energy. For IC, seed photons include (1) 
the cosmic microwave background (CMB), (2) the far-infrared 
dust emission (FIR) with the temperature of 30 K and energy 
density of 0.5eV cm ^, and (3) the near-infrared stellar 
emission (NIR) with the temperature of 5000 K and energy 
density of l.0eV cm ^. The best-fit results are shown in 
Figure 5 and in Table 2. An electron injection spectral index of 
Qe = 2.01 is required to produce the gamma-ray emission 
above 100 MeV. The total energy of the injection electrons is 
about 3.08 x 10° erg. 


^ https: / /naima.readthedocs.io /en /latest /index.html 
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Figure 4. Spectral energy distribution of SNR G332.5-5.6. Data points are plotted with 95% errors when TS > 5, while the upper limits are plotted at the 95% 
confidence levels when TS < 5. Black solid line represents the best fit while the black dotted lines are the errors at 95% confidence levels. 
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Figure 5. Theoretical expectation of the leptonic model for GeV SED of SNR G332.5-5.6. Total IC components contribute to gamma-ray emission detected by Fermi- 


LAT. Here, the electron injection spectral index is a, = 2.01. 


3.2. Hadronic Model interactions followed by pion decay (PD). In addition to 
gamma-rays, secondary electrons and positrons are produced in 
the PD process, which results in an inverse Compton emission 
(Secondary IC). Thus, the PD component plus a secondary IC 


We considered a one-zone hadronic model (Kamae et al. 
2006; Tang et al. 2017; Mo et al. 2023), in which the main 
gamma-ray production for relativistic protons are p-p 
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Figure 6. Hadronic-model fitting for GeV SED of SNR G332.5-5.6. PD and secondary IC components contribute to gamma-ray emission detected by Fermi-LAT. 


Here the proton spectral index is a, = 2.03 and the gas density is nga, = 1.00 em ?. 


component contributes for the Fermi-LAT observations. The 
distribution of the ejected protons is also assumed to follow a 
PLEC function, 


E E 
ER i A(Z) =E (- | i 


where A is the normalization, Eo is the pivot energy fixed at 
10 TeV, a, is the spectral index of the injected protons, and 
Ecutoff represents the cutoff energy. The best-fit results are 
shown in Figure 6 and in Table 2. We performed the fit and 
derived the best-fit parameters: a, = 2.03, Ecutorg = 29.82 TeV 
and the gas density of ngas = 1.00 cm~’. The total energy of the 
injection protons is about 1.32 x 10°° erg. 

In summary, both a leptonic model and a hadronic model can 
fit GeV SED of SNR G332.5-5.6 well, which share a flat 
spectrum for the ejected electrons or protons, e.g., spectral 
indices are close to —2.00. 


4. Summary and Conclusion 


In this work, we present a GeV morphology analysis of SNR 
G332.5-5.6 employing the Fermi/LAT observations from 2008 
August to 2022 December. We found that a Gaussian-disk 
spatial template with a radius of 1°06 is the preferred one with 
respect to others, e.g., a point source, a uniform disk and an 
extended template from the radio observations. Employing the 
new Gaussian disk, the GeV lightcurve for SNR G332.5-5.6 


shows a significant variability. The gamma-ray observations 
could be explained well either by a leptonic model or a 
hadronic model, in which the flat spectrum for the ejected 
electrons or protons is required. More future observations in 
other bands could shed light on the nature of broadband 
emission from SNR G332.5-5.6. 
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Appendix 
Radio Observations and the GLEAM Spatial 
Template 


The GLEAM spatial template is derived from radio 
observations of Galactic and Extragalactic Allsky MWA 
Survey between 170 and 231 Mhz(Wayth et al. 2015), the 
Gaussian smoothing is performed with bin size of 0°01, which 
can be found in Figure Al. 
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Figure A1. GLEAM spatial template (1? x 1?) from the ROI center, which is derived from radio observations of Galactic and Extragalactic Allsky MWA Survey 
between 170 and 231 Mhz (Wayth et al. 2015). The Gaussian smoothing is used with bin size of 0°01. 
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